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Abstract 

We use cosmic microwave background and large scale structure data to test a broad and 
physically well-motivated class of inflationary models: those with flat tree-level potentials 
(typical in supersymmetry). The non-trivial features of the potential arise from radiative 
corrections which give a simple logarithmic dependence on the inflaton fleld, making the 
models very predictive. We also consider a modifled scenario with new physics beyond 
a certain high-energy cut-off showing up as non-renormalizable operators (NRO) in the 
inflaton fleld. We flnd that both kinds of models flt remarkably well CMB and LSS data, 
with very few free parameters. Besides, a large part of these models naturally predict a 
reasonable number of e-folds. A robust feature of these scenarios is the smallness of tensor 
perturbations (r < 10"'^). The NRO case can give a sizeable running of the spectral index 
while achieving a sufficient number of e-folds. We use Bayesian model comparison tools 
to assess the relative performance of the models. We believe that these scenarios can be 
considered as a standard physical class of inflationary models, on a similar footing with 
monomial potentials. 
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1 Introduction 



The WMAP 3-years data [1,2] on the cosmic microwave background radiation (CMB) rep- 
resent a milestone in the history of observational cosmology due to their great precision. 
The challenge for theoretical cosmology and particle physics is to understand these data 
within some theoretical framework and to make predictions to be tested by future obser- 
vations. In this respect, inflation stands as the most successful and promising theoretical 
scenario. Proposed in 1981 as a solution to the flatness and horizon problems [3] it was 
later shown to predict an almost scale invariant spectrum of matter perturbations [4-8], 
making it a very successful idea. However, we still lack a concrete description of the 
mechanism producing inflation in the very early universe. The simplest implementation 
postulates a real scalar field, 0, whose homogeneous part drives inflation during its way 
towards the minimum of its potential, V{(j)). Therefore, to provide a description of the 
inflationary epoch of the universe we need to elucidate the functional form of the inflaton 
potential [9]. In this respect it is fair to say that, despite interesting developments, there 
is still no convincing inflationary scheme based on particle physics. 

The main goal of this work is to compare a broad class of models, which are physi- 
cally well motivated, with CMB and large scale structure (LSS) data. As shown below, 
these models predict a scalar spectral index with non-constant slope, in contrast with the 
ansatze used in standard analyses. Interestingly, these models are very predictive: the 
number of independent parameters is comparable to the free parameters used in more 
phenomenological approaches. 

This paper is organized as follows. In section 2 we review the standard way of com- 
paring inflationary potentials with cosmological data through parametrizations of the 
primordial scalar and tensor power spectra that involve the spectral indices and their 
derivatives (up to a given order). In section 3 we introduce the class of models we are 
interested in, namely, models with flat tree-level potentials which are generically lifted in 
a controlled way by radiative corrections. This type of models is then extended by taking 
into account possible effects from non-renormalizable operators that appear after integrat- 
ing out some heavy new physics. In both cases we discuss the physical motivation and 
significance of such models. Section 4 is devoted to the derivation of the power spectra for 
the two classes of models as a function of suitable model parameters. The data analysis 
procedure is discussed in section 5; we present parameter constraints within each class 
of models from our analysis in section [B] while section 7 is devoted to the issue of model 
comparison, assessing which model is in better agreement with the data. Finally, our 
conclusions are presented in section 8. In addition, the appendix contains some technical 
details relevant to the analysis of the model with non-renormalizable operators. 



2 Comparing inflationary potentials with data 

The usual strategy to probe inflationary potentials is the comparison of theoretical pre- 
dictions for the spectrum of primordial (scalar and tensor) cosmological perturbations 
[Ps{k) and Pt{k) respectively, where k stands for the space wave-number] with CMB and 
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LSS data. This is usually done in the slow-roll approximation [10]. Generically, in this 
situation |e, ?7,^| <C 1, with 



e 



where primes denote derivatives with respect to the inflaton (j) and Mp = nip/ v^Stt (being 
rrip ~ 10^^ GeV) is the reduced Planck mass. In this approximation 

3 V 

PJk) ~ • (2.3) 

Thus the tensor to scalar ratio simply reads 

p 

r = — ^lQe. (2.4) 

The relationship between the inverse distance scale k and the inflaton fleld (p is given in 
slow-roll by 

-^.-M,v/2i. (2.5) 

The number of e- folds between the beginning of inflation {t^) and the end (tg) is 

I ft" V 1 Z"'^* 1 

where we have used the slow-roll approximation to change variables from t to 0. The 
usual requirement to solve the horizon problem is Ai'e — 50 — 60 although the precise 
number depends on the details of the post-inflationary epoch [11]. Of course, the total 
number of e-folds may be much larger, but only 50-60 e-folds are needed and actually 
WMAP only probes directly the flrst ~ 7 e-folds, which can be extended to ~ 10 by using 
LSS data and other observations. Thus the subscript above denotes the time at 50-60 
e-folds before the end of inflation. On the other hand, there might be several episodes of 
inflation, in order to achieve the required 50-60 e-folds. In that case, WMAP would be 
only sensitive to the first (s) of these episodes after t^,. 

An important quantity when performing fits to the observations is the scalar spectral 
index, n, which describes the variation of Ps with /c, 

c/ In P, , , 

n-l = — — ^ ~ 2r/ - 6e . (2.7) 
am A; 

The spectral index itself may change with k: 



(2.8) 



The right-hand sides of (12. 7p and (12.81) correspond to the slow-roll approximation. Similar 
equations, involving the tensor spectral index nt{k), can be written for the variation of 
Pt{k) with k. Summarizing, observational information about Pg (and thus about n and 
dn/dlnk) gives information about V and its derivatives, V, V", V" and so on, providing 
a link between theory and experiment. 

Standard Parametrizations of the Power Spectrum 

A common approach to determining the power spectrum from CMB and LSS data is 
to perform a Taylor expansion of XTvPgik:) and lnPf(/c) in In/c/fco around zero, where k^ is 
a pivotal scale (from now on fixed to k^ = 0.002 Mpc~ ), 

k 1 dn / I- \ 2 



\nPs{k) = \nPs{ko) + [n{ko)-l]\n— + 
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fco 



InA) (2.9) 



\nPt{k) = \nPtiko)+nt{ko)\n^ + -- - . (2.10) 

ko 

Note from equations (12. 2p and (12.41) it follows that, at first order in slow-roll [12]: 

r{ko) = ^ ^ -8n,{ko) . (2.11) 

Usually, one does not go beyond the order shown in the equations (12.91) and (I2.10p . In par- 
ticular, a running of the tensor spectral index has not been considered because presently 
the tensor contribution to the spectrum is only weakly constrained by the data. Therefore 
one typically fits four independent parameters, namely: {In Pg, n,dn/d In k, r}ko and often 
the running of the spectral index, dn/dlnk\kQ, and the tensor to scalar ratio, r{ko), are 
set to zero. The improvement in the fit obtained when introducing the latter parameters 
is neither large enough to be considered as a strong indication for their presence nor small 
enough to be considered as irrelevant. In general, the issue of whether an extra parameter 
is needed or not is a difficult one and has to be addressed with care, see [13, 14] for a 
discussion. We return to this question in Section [3 

The Taylor expansion above has been used by the WMAP collaboration. Assuming a 
ACDM universe and setting dn/d\nk = r = 0, WMAP 3-yr data alone give [15] 

n = 0.958 ±0.016 , (68% c.l.) (2.12) 

Including a constant dn/dlnk the result of the analysis is 

n{ko) = 1.050^^58, (68% c.l.) (2.13) 
^ = -0.055+°;°^; (68%c.l.). (2.14) 

The absolute magnitude of Pg depends slightly on the inflationary model, but roughly one 
flnds Ps{ko) ~ 2 X 10~^. The WMAP collaboration used these results to probe monomial 
potentials V{(f)) oc 0° (with a = 2,4) [2]. These models predict a negligible running of 
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n, so they are well approximated by equations (12.91) and fl2.10p with dn/dlnk = 0. The 
fits seem to exclude a = A [16] and any other higher monomial power. The quadratic 
case, V = |m^0^, works quite well, although it requires very large values of the inflaton 
field, ~ 14 Mp. For the purposes of our discussion it must be stressed that the simple 
functional forms assumed previously for Ps{k) and Pt{k) may not be accurate enough to 
describe the actual power spectrum of other inflaton potentials which are well motivated 
physically. In this sense, although (12. 9p and (I2.10p can be useful as phenomenological 
approximations, it is important to be open to other parametrizations. In this paper we 
will show explicit examples of this. At the end of the day, the best fit together with the 
best physical motivation will determine the preferred functional form. 



3 A broad class of models: Flat Tree-Level Potentials 

We will consider models that have "/Zai tree-level potentials" 0, i.e. 

Vtree{4>) = Ptrco = COUStaut . (3.1) 

Then, the potential derivatives V, V",... arise from the radiative corrections to V. These 
potentials appear typically in supersymmetric (SUSY) theories: V^trec^'*' ordinarily has 
plenty of accidental fiat directions. A familiar example of this is the minimal supersym- 
metric standard model (MSSM). Such accidental flatness is broken by radiative corrections 
since there is no symmetry protecting it. Generically, at one- loop 

m=P + /5 1n^, (3.2) 

where Q is the renormalization scale (which might have absorbed finite pieces) and m(0) 
is the most relevant 0-dependent mass in the spectrum. Note that p depends implicitly 
on Q through its renormalization group equation (RGE) and that the Q-invariance of the 
effective potential implies 

at one-loop. From now on we will assume /3 > 0, which is the usual situation, though the 
opposite case is also possible and the analysis is similar. 

Since these "almost fiat" directions are so common in SUSY scenarios, they are natural 
candidates to drive inflation, provided the potential stores large enough energy density. As 
a matter of fact, particular examples of these approximate flat directions have been used 
in the literature to implement inflation, e.g. the first D-hybrid inflation model belongs to 
this class [18]. The important point is that, whatever the model considered, the slope of 
V{(j)), and thus the dynamics of inflation, are determined by radiative corrections. Since 
the latter have a very generic functional form (logarithmic), it is possible to make very 



For a more complete diseussion of the theoretical aspects of these models see [17] 
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model-independent predictions without relying on a particular model [17]. Next we work 
out these statements in a more detailed way. 

The leading-log approximation (which amounts to summing up the leading-log con- 
tributions to all loops) is implemented in this context by simply taking Q = m{(j)). This 
choice eliminates the potentially large (and thus dangerous) logarithms, improving the 
convergence of the perturbative expansion. Then 

V{<P) ^ p[Q = m(0)] . (3.4) 

In general, one expects m^(0) = + c^^^, where M does not depend on 0, and c is 
some coupling constant (which depends on Q according to its own RGE). Normally one 
considers a range of 0— values where either the constant piece (if it exists) or the 
0— dependent part dominate. We assume we are in the second case. Hence, we will ignore 
the possible presence of M and take Q = m{(f)) = ccp. 



Logarithmic regime (LOG) 



In the regime of very small coupling constants one has d(3/dlnQ ^ f3 since the former is 
higher order in the couplings and has a loop suppression factor. Then we can consider f] 
as constant in the range of Q cx of interest (which is never too wide). Now the scalar 
potential 03.41) can be easily written in terms of its value at 0o, the value of the inflaton 
at fco, integrating equation 03.31) : 

V{4>)=po + f3 In-^, (3.5) 

00 

where po = p{(po). Note that 03.51) can also be obtained from 03.21) by choosing Q = c0o. 

Since ln</) = limQ,_»o a~^(0" — 1) the logarithmic potential 03.51) can be considered in 
many respects as a monomial potential with a = 0. In particular all the derivatives, 
which are crucial for the cosmology of the model, follow that pattern. On the other hand, 
as argued above, this potential is physically as well motivated as the monomial forms. 



Logarithmic regime -|- non-renormalizable operator (LOG-|-NRO) 

If there is a scale of new physics, M, higher than the scales relevant to inflation (i.e. 

S> 0^), the new physics will generically show up in the effective theory at lower scales 
as non-renormalizable operators (NROs) of the light fields, suppressed by inverse powers 
of M. Due to the suppression factor, the impact of the NROs in the physics at low scales is 
usually very small. However, if the NRO has characteristics not shared by the low-energy 
physics, its effect may be significant (as happens with higher-dimension operators that 
mediate proton decay or give a Majorana mass to neutrinos). In our case, the new physics 
does not need to respect the accidental fiat directions of the effective theory. Thus one 
expects the inflaton potential 03.51) to become 

V{<p) = po + (3\nj- + <jy'^. (3.6) 



The first two terms correspond to the generic one-loop potential in the small-couphng 
regime while the last term is a non-renormalizable operator left in the low-energy theory 
after integrating out some unspecified physics at the high scale M. Notice that this scale 
absorbs any possible coupling in front of the operator. Of course V{(j)) may contain other 
NROs of different order. Here we assume that the one shown in equation fl3.6p is the 
lowest order one, and thus the dominant one. The sign and power we have assumed for 
this NRO guarantee the stability of the potential. Notice also that an even power for this 
operator is what one expects generically in supersymmetric theories. An explicit example 
of this is given in [17]. Apart from that, the potential ( 13. 6p is completely general and 
therefore the analysis is quite model-independent. 



4 Primordial spectra in the slow-roll approximation 
4.1 Logarithmic regime (LOG) 

From ^(0), as given in equation (13. 5p . the first three slow- roll parameters (12.11) in the 
LOG scenario are 

1 e 

^-o^'^' ^--?-f--2-, e^2r/2, (4.1) 

2 (p^ (p^ q 

where we have introduced the quantity 

q^-- (4.2) 

Po 

Note that g > since we have assumed /5 > 0. Satisfying the slow- roll condition e ^ 1, 
requires g <C 1 for < Mp. This implies in turn e <^ |?7|. We give a more quantitative 
evaluation of this hierarchy in Section 5. 

The number of e-folds between time to (i.e. the time when (p = (po) and the end of 
inflation, tf., can be easily computed using e, as given by equation (14. ip . in the usual 
expression (12.61) : 



1 



[r]{(po) Vi<l>e)\ 2r/(0o) 2qM^ 



^ ^ '^o - iV" , (4.3) 



where we have used the fact that inflation comes to an end when t] ~ 0{1) (it could 
end before that time if inflation is interrupted by other mechanisms, like a waterfall 
field in hybrid models). The parameter defined above, besides giving an excellent 
approximation to A'^p(to t^), will play a relevant role when performing the fits to the 
data. 



Now the (p — k connection, equation (12.50 . can be integrated at first order in q, giving 

1 k 



ii(i-4'4i. (") 



e 
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where we have used (14.31) . Note that increasing corresponds to decreasing k so that 
the scales probed by WMAP correspond to the highest values of during its slow-roll 
towards the origin. Now we can straightforwardly evaluate Ps{k) from equation (12. 2p . For 
the purpose of comparing the model with the data, it is convenient to write Pg in terms 
of P° = Ps{ko) using the general expression 

In = In P° + 3 In ^ - 2 In ^ . (4.5) 

Using equations (13.51) and (14.41) and expanding in q we find, at first order, 

lnP,(A:) = lnPi'+(^l + ^g)ln(^l-^lnA^ , (4.6) 



where, from equation (12.21) . 



pO = _± ^0*^0 (A 7) 

' 127r2 M6/?2 ■ ^ • ^ 

The same result can be obtained by integrating the slow- roll equation (12. 7p . 

Similarly, the spectrum of tensor perturbations, Pt^k), can be obtained from equa- 
tion (O 



At the same level of approximation, the tensor-to-scalar ratio (12. 4p reads 



(4.8) 



Let us now count the number of independent parameters. The power spectra, Ps{k) 
and Pt{k), contain three independent parameters, {Pf, q, N^}, which are combinations 
of the initial parameters {(pQ, po, (3}. Incidentally note that the scalar potential (13. 5p is a 
function of just two combinations of parameters, but a third one appears in the conversion 
of (j) into k through equation (14.40 . Actually, the q term in (14.60 is subdominant because 
g -C 1. Removing it from the expression is a good approximation and eliminates one 
parameter. So the expression of Ps{k) contains basically two parameters. This is to be 
compared with the three parameters (two if the running of n is set to zero) of the simple 
standard parametrization (12. 9p . As a consequence this LOG scenario is highly predictive. 
On the other hand, the fits to WMAP data prefer Pt{k) <^ Ps{k), which means that Pj 
turns out to be scarcely important in the fit, and so the number of relevant parameters 
continues to be two. Indeed, from equation (14.90 and g <^ 1 we do expect by construction 
Ptik) <^ Ps{k), something that cannot be postulated from the simple parametrizations 
(12. 9p and (I2.10p . which contain the additional parameter Ptiko) or, equivalently, r (unless 
r is set to zero by hand). 

Equations (14.60 and (14.80 summarize the predictions of models with flat tree-level 
potentials in the small coupling regime. It is possible to gain some intuition about the 
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performance of these functional forms as follows. As mentioned above, e <^ |?7|, which 
means n — 1 ~ 2//. Therefore [see equation (14. ip ]. 

-2e = -(r^-l)^ (4.10) 



dink 



As a consequence, dn/dlnk is negative, as suggested by the data, though its value tends 
to be quite small. In fact the sign of n — 1 cannot change along the inflationary process 
in this class of models. Since the dependence of n on the scale is weak, and therefore 
n ^ constant, we can expect a fit similar to the one obtained by using the standard 
parametrization of equation (12. 9p with dn/dlnk = 0, leading to no ~ 0.95. We will see in 
Section 5 that this is indeed the case. 

Equation (I4.10p can be written in an integrated form as 

" = ' - ATO - ln(A:/A:o) ' ^''''^ 

where we have used equation (14.30 . Equation (14. lip is the prediction for the spectral 
index in scenarios with fiat tree-level potential in the regime of small coupling. It can be 
compared with the n = constant or dn/dlnk = constant assumptions made in standard 
analyses. Note that is the only independent parameter in (14. lip and has a precise 
physical meaning. 



4.2 Logarithmic regime and non— renormalizable operator 
(LOG + NRO) 

In this scenario the derivatives of the potential V{(f)) [equation (13. 6p ] with respect to 
read 



V"{<P) 
V"'{<P) 



6 6"^^ 



^ + 2(A^ + 2)(2A^ + 3)02 



l2N 



M2^ 



>3 



+ 4(A^ + 2)(2A^ + 3)(A^ 



l2N 



M2N 



The corresponding expressions for e, r] are 



1M2 r 



V 



p_ 

2 02 

M2 

p 



/^4 a2N 



-£ + 2(iV + 2)(2/V + 3)r-^ 



/,4 jJ2N 



(4.12) 



(4.13) 



In consequence the NRO can have a significant impact on inflation when the small number 
(0/M)2^ is comparable in size to /3/0^. It is also immediate to realize from equations 
(I4.12[ I4.13P that, for sufficiently large 0, the higher derivatives V" , V" (and thus r/, ^) can 
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receive a large contribution from the NRO while the contribution to V (and thus e) is 
much less significant, thanks to the additional (2A^ + 3) and (2A^ + 3)(2A^ + 2) factors in 
V", V". Therefore, one expects modifications of the spectral index, n, and its running 
dn/dlnk (which depend on r],^), especially at the initial scales (very small k and thus 
large 0), and much smaller changes in e (and hence in A'^g)- 

Let us now calculate the expression of the power spectrum Ps{k) in this scenario. As 
in the previous subsection, we start with the general expression (14.51) where V{(f)) and 
V'{(l)) are given now by equations (13. 6p and (I4.12p respectively, and 



127r2M2[\/'(0o)]^ 



(4.14) 



Next we have to convert </> into k by integrating equation (12.51) . This is done in the 
Appendix in an exact way. The corresponding formula, equation ( lA.ll) . is rather cumber- 
some. However, using |g| ^ 1, one can write the much simpler but extremely accurate 
expression 



1 ^ 

ko 



1 



A^- 



1; 



iVH-3_ 
A^ + 2' 



N+2 



^=1 



(4.15) 



where 2-^1 (a, b; c; z) is the Gauss' hypergeometric function [19], ip is just a dummy variable 
and 



2(A^ + 2)- 



(4.16) 



/5 ■ 

Note that N° is still defined as in equation (14.31) and gives a good estimate of the number 
of e-folds A^e(^o ~^ te)- It is possible to invert (14.151) numerically to get </> = </'(/c). Plugging 
4>{k) into equation (14.51) we obtain Ps{k). This is the procedure we have followed in doing 
the fits. Using 0(/c) we can also obtain other quantities of interest as functions of k, e.g. 
the spectral index n ^ 1 + 2?7 — 6e or the tensor to scalar ratio r ~ 16e. 

In order to get some intuition about the shapes of Ps{k) and n{k) it is convenient 
to derive an analytical approximation to the previous numerical procedure. Actually the 
numerical part just comes from the (pto k conversion, i.e. the integration of equation (12.51) . 
This equation depends on e (on V and V) but not on higher derivatives of V, which are 
the ones most affected by the presence of the NRO, as discussed after equation (14.121) . 
Therefore it is sensible to use here the value of e when the NRO is switched off, i.e. that 
of equation (14. ip . Then the (p to k relation is still given by equation (14. 4p . Substituting 
this in the general expression for Ps{(f)), equation (14.50 . and expanding at first order in 
the NRO contributions, one get^ 



In 



Psjk) 
P? 



-g ) In I 1 



1 k 



N + 2 



1 k 



N+2 



(4.17) 



This approximate formula gives Ps{k) with a maximum error of ^ 13% in the most extreme cases 
although typically is much better. Anyway, wc remark that in the fit we evaluate Ps{k) numerically using 
equation (|4.15p . 
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where 



7 = <^ 2iV + 3 



2{N + 2) 



-3 



o ^ aN+2 aN+2 



The approximate equahty in the last expression is justified by the smallness of q. Alter- 
natively, one can expand n — 1 ~ 2?7 — 6e at first order in the NRO and approximate again 
the (p to k conversion by equation fl4.4p . One obtains 



3 \ 1 / 1 . k\-' /. 1 , *-^"+' 



„W-l.-^l + -,j-^^l-^l„-j . (4.19) 

Then the direct integration of ( 14.19^ gives back expression (]4.17p . Equation ( 14.19^ cor- 
responds to a running n{k) with non-constant slope, departing from the assumption of 
analyses done using the standard parametrization. Unlike in the LOG scenario, in this 
case the running is not constrained to be very small. 

It is also worth mentioning that, due to the positivity of 1 — (l/N^) Ink/ko, the sign of 
the LOG and the NRO contributions to {ra(A;) — 1, dn/dlnk, (Pn/dln^ k} are { — , — , — } 
and {+, — , +} respectively [see equation fl4.19l) ]. Since a sizeable running at low k requires 
a dominant NRO contribution, we can conclude from equation fl4.19p that in that case 
the sign of the second derivative will be positive, although for large enough k it will turn 
to negative as the LOG part becomes dominant. 

In a similar way one can obtain expressions for Pt{k) [or equivalently r{k)] starting with 
the general equations (12.30 or ( 12.40 . In particular, the previous analytical approximation 
gives in this case 

4g / 1 , k\-^ 8gA^+V q \f 1 , k^^^ , , 



Let us finally count the number of independent parameters. From expressions (14.171) 
and (I4.20p we see that the spectrum of primordial perturbations depends upon five param- 
eters {-Pf, Ai'g , g, A, iV}, which are combinations of the five parameters {po, 0O) ^! ^} 
appearing in the scalar potential (13.61) . Thus the LOG+NRO has two more parameters 
than the LOG model. Again, as in the LOG case, the smallness of q implies that Ps(A;) is 
nearly independent of q and, besides, the tensor spectrum is much less important than its 
scalar counterpart. Hence, q will be irrelevant for a broad range of values in the fit to the 
data. In practice the primordial spectrum depends essentially on four parameters which 
become just three if we consider iV to be a fixed integer. Again, this is to be compared 
with the four parameters of the simplest standard parametrizations (12.90 and (12.101) . In 
consequence, the LOG-I-NRO scenario is still highly predictive. 

To illustrate the shapes of the scalar power spectrum in the different scenarios dis- 
cussed in this section we plot Ps(/c) from (I4.15P in Figure [H The parameters for the 
models are chosen to be the best-fit values given in Section 6 and discussed later in the 
text. 
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Figure 1: Primordial power spectra in the standard parametrization (dashed hnes with 
Taylor expansion up to second order, i.e. running of running) and as predicted by the 
LOG and LOG+NRO {N = 2) scenarios. Parameters for each case are the best-fit ones 
(given in Tables [1], [2] and [3]) . 
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5 Data analysis procedure 



In order to constrain the parameters of the two scenarios introduced above, by comparing 
their predictions with CMB and LSS data, we use a modified version of the cosmomc 
package [20] with a suitable parametrization of the expressions for Ps{k) and Pt{k). In the 
next two subsections we describe the detailed procedure for each case separately. We ran 
six chains for each model, gathering 3x10^ samples per chain, using the default Metropolis 
algorithm to sample the parameter space. We discard a burn-in period encompassing 
the first 10'^ samples in each chain and we employ the usual Raftery and Lewis mixing 
criterium [21], requiring i? — 1 < 0.1 for the merged chains. 

We consider a fiat cosmology, taking fiat priors on the cosmological parameters flbh'^, 
VLf.h? (the baryon and CDM density, respectively), 9* (the ratio of the distance to the 
last scattering surface to the sound horizon) and r (the optical depth to reionization) . 
We assume dark energy in the form of a cosmological constant [w = —1) and 3 species of 
massless neutrinos. The parametrizations of the primordial power spectrum are discussed 
in detail below. 

We use a combination of CMB and LSS data to constrain the parameters of the models. 
The main reason for doing so is that it helps to break parameter degeneracies [22]. In 
particular, we use the latest WMAP [23], ACBAR [24], CBI [25] and BOOMERANG [26] 
CMB data set releases and the results from the SDSS galaxy survey [27]. We also add the 
Hubble Space Telescope measurement of the Hubble constant Hq = 72±8 km/s/Mpc [28]. 



5.1 Logarithmic regime (LOG) 



As discussed in Section 14.11 the LOG scenario is described by the three independent 
parameters {p, ^o}, appearing in the scalar potential (13.51) . However, for the purpose 
of comparing the model with data, it is more appropriate to use the following set: 

F^OG = UnP^,N'^,q}, (5.1) 

which are related to the potential parameters by the relations (14. 2p . (14. 3 p and (14. 7p . The 
inverse transformations are given by 



= ^M' , (5.2) 



The reasons for preferring the set (15. ip over the original potential parameters are the 
following. First, the fit to WMAP data is very sensitive to the value of InP^ at the pivotal 
scale, which makes it very convenient to use lnP° as one of the parameters. Second, 
appears explicitly in the expressions of Ps{k) and Pt{k) [see equations (14. 60 and (14. Sp . 
In addition, A'^^ has a clear physical meaning, since for small q it simply expresses the 
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number of e-folds since the time when ko crosses out the horizon until the end of inflation. 
This also makes it possible to impose on it a physically motivated prior for the number 
of e-folds, as required to solve the homogeneity and flatness problems. Furthermore, we 
note from equation fl4.1ip that A'^^ and the spectral index at ko are simply related. 



Finally, the third parameter, g, does also appear explicitly in the expressions for the 
spectra. As argued in subsection 14.11 we expect g <C 1, implying that the scalar pri- 
mordial spectrum depends essentially only on {In Pg,N^}, while the tensor spectrum 
Pt{k) ^ 16ePs{k) is suppressed (and much less important for the fit). Therefore it is con- 
venient to choose g as a parameter for the fit in order to single this effect out. From 
the above discussion. In Pg,N^ will be well-determined by the observable properties of 
the power spectrum, and therefore it is appropriate to impose fiat priors on them, which 
corresponds to the assumption that they are location parameters. 

The relations between Plog and the potential parameters is non-linear and so one 
expects volume effects coming from the Jacobian of the transformation that will in general 
make the marginalized constraints on the potential parameters sensitive to the choice 
of priors. Furthermore, as argued above, only two combinations of parameters of the 
potential are going to be well determined by the data. The constraints on these "principal 
directions" in the potential parameter space are however essentially prior-independent, 
as we discuss in detail in Section 16.21 

Let us now focus on the physical constraints on the parameter space spanned by Plog- 
The evolution equations of the classical value of the field are based on General Relativity. 
To prevent effects of quantum gravity from becoming important, we conservatively require 
the energy density to satisfy 



Similarly, it is sensible to keep the infiaton field below the Planck scale. Note in particular 
that, at least in this framework, the renormalization scale Q is to be identified with the 
value of the infiaton, in order to maintain the radiative corrections under control, and 
obviously the RGE are only reliable for Q below the Planck scale. Thus we also require, 
conservatively. 



Notice that, since the infiaton rolls towards zero, if the above condition is satisfied for 
00 it will automatically be satisfied for smaller values of (j), as well. For larger values of 
0, imposing equation (15. 5p easily guarantees that they are well below nip, since there are 
very few e-folds before ko, and they correspond to a short range of 0-values. 

Moreover, one must ensure that the slow-roll approximation is fulfilled, which means 
that we require 




(5.4) 



00 < Mj 



(5.5) 



e < 1, 
r]\ < 1. 



(5.6) 
(5.7) 



13 



For simplicity we impose these conditions at ko and this automatically ensures that the 
slow-roll is not violated for smaller values of k, which means greater values of 0. Therefore 
the slow-roll will be guaranteed at the scale kobs = 10~^Mpc~^, which is roughly the size 
of the observable universe. On the other hand, larger values of k are probed and the 
slow-roll parameters grow as goes to zero. Taking into account that the largest relevant 
multipole is about 3000, one gets a maximum k around k^ax = 0.1 Mpc~^ [29]. Using 
equation (14.41) . we have checked that slow- roll for such large value of k is indeed satisfied 
by the samples in our Markov chains. The slow-roll condition on r] is equivalent to 

2iV° > 1 , (5.8) 

while the one on e leads to 

AN^ > q. (5.9) 
On the other hand, the inequality (15. 5p implies 

2gA^°<l, (5.10) 

which together with equation (15. 8p implies 

g<l, (5.11) 

as anticipated in Subsection 14.11 

We found that samples in the Markov chains that fulfill the condition (15.101) . automat- 
ically satisfy also conditions (15.81) . (15. lip and (15.40 . This can be understood as follows. 
As discussed in Subsection 4.1, we expect a value for uq ~ 0.95, similar to the simplest 
fit (I2.12p . Then equation (15. 3p implies (15.80 . Moreover, the value of n, coupled with the 
physical prior (15. 5p . translates into an upper bound on q: since n — 1 ^ —2rj ~ 2g(Mp/0o)^5 
for 00 < one gets q < (1— n)/2, and therefore condition (15. lip holds. Incidentally, this 
upper bound on q implies that the contribution of the tensor part of the power spectrum 
must be necessarily small: r = 16e < 2(1 — n)^. Finally, (15.41) is granted by the smallness 
of P°. Notice also that the condition (15. 6p on e is, in practice, irrelevant because (15.111) 
ensures that e < |?7|, as can also be read off directly from (14. ip . In consequence, condition 
(I5.10p remains the only non-trivial constraint. 

In summary, we take Plog [equation (15. ip ] as the set of independent parameters, 
imposing fiat priors on them and enforcing the constraint (I5.10p . We then compute the 
scalar and tensor contributions to the primordial spectrum via the expressions (14. 6 p and 

iH). 

As mentioned above, one of the reasons for choosing as an independent parameter 
is its direct physical interpretation as the number of e-folds. In fact, we have a strong 
theoretical prejudice about its value, which should 50. We have taken into account 
this fact by performing two different analyses. The first one imposes a flat prior on A'^^, 



^Recall however that its value could be less than 50 if there are subsequent episodes of inflation [see 
discussion after equation (|2.6p ]. On the other hand, the parameter A^" [defined in eq. (|4.3[) ] could be 
larger than 50 if inflation is interrupted e.g. by a waterfall condition in hybrid models [see discussion 
after equation (|4.3p ]. 
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therefore assuming no prejudice about its value and leaving the data to constrain it. In 
the second case we enforce the theoretical requirement by imposing a Gaussian prior on 
A'^^ centered on 50 with a standard deviation of 5. The details of these two fits and the 
results are discussed below, in Subsection 16.21 

5.2 Logarithmic regime and non— renormalizable operator 
(LOG + NRO) 

Concerning the LOG+NRO scenario, for practical reasons it is convenient to work with 
the set of independent parameters 

Plog+nro = { In P°, iVO, g , A , iV } , (5.12) 

where A was defined in equation (14.161) . instead of the parameters {p, /?, 0O) M, N} of the 
scalar potential (13.51) . The relationships between Plog+nro and the potential parameters 
are given in appendix El 

The convenience and significance of the first three parameters in Plog+nro are the 
same as in the LOG scenario. However, the interpretation of A''^ as the number of e-folds 
between ko and the end of inflation is now less accurate since there are NRO corrections, 
although it is still a good approximation. This is also true for the connection between 
and the spectral index no: the relation (15. 3p becomes now 

/ 3 \ 1 1 

This expression tells us that for not too large we should expect A to be bounded 
from above by some number close to unity because otherwise no can become substantially 
different from 1 (especially for high values of A^), thus violating the slow-roll conditions. 
This is also illustrated in Figure [2] which shows contour plots for the scalar spectral index 
and its running at at lowest order in slow-roll as functions of A and A'^^ for N = 2 
and A^ = 10. It is worth remarking here that since we are dealing with scale dependent 
quantities the appearance of these graphs would change if we made them at a different k. 
Figure 2 allows to realize that, in the context of the LOG-(-NRO scenario, it is possible 
to have simultaneously a sizeable running and a reasonable number of e-folds. 

Concerning the physical limits in parameter space, we must take into account the 
presence of the scale of new physics M. The role played by the Planck mass on the LOG 
scenario corresponds now to M. To keep the validity of the effective potential (13.61) the 
infiaton must evolve well below that scale, which should be itself smaller than the Planck 
scale. Thus, we impose the following conservative limits: 

2^obs <M<Mp, (5.14) 

p<M\ (5.15) 

Notice that we set the first constraint at (pobs = <Pikobs) in order to ensure it for any value 
of in the observable range. 
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Figure 2: Values of the scalar spectral index (red lines) and its running (blue dashed 
lines) at /cq computed at lowest order in slow roll for the LOG+NRO class of models with 
N as indicated. 



As in the LOG scenario, |?7| ^ e, so ?7 is the relevant parameter for the breakdown of 
the slow-roll. However, unlike in the LOG case, due to the NRO the absolute value of i] 
grows with sufficiently large 0. Therefore we must ensure the fulfilment of slow-roll not 
only at kmax but also at kohs- This guarantees that any point in between will satisfy the 
slow-roll conditions as well. So, we reject in the Monte Carlo process those points such 
that 

|?7(fco6s)| > ?7lim , (5.16) 

or 

Vikmax) > Vlim, (5.17) 

being ?7iim a limiting value (smaller than 1) that we set at the beginning of the run. At 
the end we check that rjn^ was indeed well chosen to ensure the validity of the slow-roll 
approximation. In practice we work with ?7iim = 0.2 which is a rather conservative value. 
We have checked that the change in the results is negligible if instead we use ?7iim = 0.5. 

On the other hand, we expect theoretically that the parameter should be in the 
range from 1 to (9(10). It is very common that fiat directions in supersymmetric models 
are only lifted by NROs at very high order (as it is the case of the MSSM [30,31]). This is 
also very common in D=4 string compactifications due to stringy selection rules [32,33]. 
For further details see [17]. In Section [6] we discuss in detail two representative cases, 
which reasonably encompass the range of values for A^ (A^ = 2 and A^ = 10), and we 
comment on the qualitative behaviour for values of A^ in-between. 

As for the LOG scenario, we can consider as a free parameter with a fiat prior on 
it or we can alternatively constrain it to be around 50. We have performed the two types 
of fit. 
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Finally, we can anticipate theoretically the appearance of some strong bounds on the 
parameters of the model when performing the fits. First note that the observed power 
spectrum normalization P^" ~ 2 x 10^^ implies through equation (12.21) the smallness of 
Po/Mp. More precisely 

(5.18) 



p 



where the subscript "0" indicates evaluation at the pivotal scale. On the other hand, 
from equations (14.131) we note that the smallness of \r]\ (to preserve the slow-roll) im- 
plies that the two contributions within the square brackets (i.e. the LOG and the NRO 
contributions) must be small separately. Otherwise one should require an unjustified fine- 
tuned cancellation between them. Actually, even with fine-tuning, one could arrange the 
parameters to produce the cancellation only at a particular (and thus k): since the 
0— dependence of the two terms is very different, at another (not too distant) value of 
the cancellation would not work, spoiling the slow- roll. Consequently, the smallness of |?7| 
implies 



I 
Po 

2N+2 



< 



< 



Vo 



1^0 I 



Po 



5 X 10"^ 



eo|'7o| 



2(A^ + 2)(2iV + 3)M4 2(A^ + 2)(2iV + 3) 



(5.19) 
(5.20) 



In the second equation we have used M < Mp and equation (15.181) . On the other hand, 
comparing the two equations (14.131) . and recalling that there cannot be fine-tuned can- 
cellations in r}, it is clear that 



Using this relation in (I5.20p we get 



(5.21) 



2N 



which, substituted into (15.191) . gives 



7U 13 



< 



5 X 10-'\rio 



4(iV + 2)(2iV + 3) 



(5.22) 



Po 



5 X 10- 



4(iV + 2)(2iV + 3)_ 



(5.23) 



This shows that q is typically small: for N = 2 [N = 10) one obtains q ^ 1.7 x 10~^ 
(g ^ 1.4 X 10^^), a conservative estimate obtained by replacing tiq 1]^^ = 0.2. In 
practice t] is substantially smaller since the bound rj < rjn^ is to be fulfilled at all k, not 
just at the pivotal scale. Similar bounds on q can be obtained using our priors together 
with the constraint M < Mp. We have checked numerically that the values obtained 
agree with the ones derived above. 
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6 Results 



6.1 Standard parametrization 

The results obtained using the standard parametrization for the primordial spectra [equa- 
tions ( 12. 9p and ( ]2.10p ] are summarized for easy reference in Table [1] both with and without 
a running spectral index [i.e. including or not the last term of equation (12. 9p ]. For later 
reference, we have also considered the next term in the Taylor expansion (12. 9p . which 
has been denoted as "running of the running" , in the last two columns of Table [H The 
table shows the best-fit parameter values, the posterior values and 68% 1-dimensional 
posterior intervals for the parameters. We also give the best-fit values for (minus twice) 
the log-likelihood, normalized with respect to the model with only a constant tilt in- 
cludeqj. Let us recall that the quantities describing the primordial spectrum are defined 
at ko = 0.002 Mpc~^. These results will be useful later on to interpret the outcomes for 
the LOG and LOG+NRO models and for the comparison with them. 

It is interesting to note that when using the standard parametrization up to second 
order a large and negative running of the running is preferred [34], which increases the 
power on large scales (see bottom panel of FigureHj), even though the fit is only marginally 
better than the case with constant running, see Table [H This is somewhat surprising: 
since in the slow-roll approximation n = 1 + 2?7 — 6e, if n departs from scale-invariance 
too quickly then rj, e or both grow up to 0{1) values, marking the end of slow-roll and the 
inflationary process at quite small k. However, in order to solve the horizon problem we 
need ~ 50 — 60 e-folds of inflation, which corresponds to the same interval in In k. This 
requires the (large and negative) dn/dhik to get suppressed at some point, suggesting 
a positive second derivative (unlike the result of the fit). This contradiction could only 
be avoided if the n{k) function changes abruptly at some point (or, maybe, if there are 
subsequent episodes of inflation). In any case, it seems clear that the current preference 
for a negative second derivative is strongly driven by the (possibly anomalous) low power 
of the large-scale multipoles. This could easily change if observations by Planck do not 
confirm the lack of large scale power observed by COBE and WMAP. However, we notice 
that from a model selection perspective even present-day data do not require a non-zero 
running of the running, as discussed in Section 7. The previous discussion will be useful 
later on to interpret the outcomes for the LOG and LOG+NRO models and for the 
comparison with them. 



6.2 Logarithmic regime (LOG) 

We denote the choice of flat prior on by LOGJF, where the top-hat distribution is 
taken in the interval 2 < A'^^ < 1000. In the second case, denoted by LOG^, we impose a 
Gaussian prior on iV°, with mean 50 and a standard deviation of 5 e-folds. Let us remind 

■''The absolute value of the log-likelihood is of little interest here and in the following. For completeness, 
we have computed the likelihood values using the WMAP3 likelihood code version v2p2 with the default 
settings regarding the offset for the log-likelihood. The best-fit value for the constant tilt model is 
-21n£ = 3614.0. 
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Table 1: 1-dimensional marginalized 68% region and best fit values results for the standard 
parametrizations (12.91) . fl2.10p with n = constant, dn/lnk = constant and d'^n/ln^k = 
constant, from left to right. 



Model 
-2 A In/: 



no runnmg 



0.0 

ID 68% Best fit 



with running 



-3.4 

ID 68% Best fit 



running of running 



-4.4 

ID 68% Best fit 



Cosmological parameters 



Qi,h^ X 10^ 
r 

i/o[Km s-^Mpc"^] 


2.23 ±0.07 2.25 
0.106 ±0.004 0.107 
1.043 ±0.003 1.042 
0.084 ±0.029 0.087 
74.3 ±2.1 73.1 


2.20l^;|^^ 2.18 
0.107 ±0.004 0.109 
1.043 ±0.003 1.043 
0.114 ±0.035 0.113 
73.1 ±2.3 72.0 


2.22 ±0.08 2.18 
0.107 ±0.004 0.107 
1.044 ±0.004 1.043 
0.106 ±0.033 0.109 
73.6 ± 2.4 72.8 


Power spectra parameters 


ln(po X 10^°) 

no 

d^n/dln^fc|^^ 


3.11 ±0.07 3.15 
0.973 ±0.019 0.961 

< 0.22 0.003 


3.00 ±0.10 3.09 
1.14ll°:°| 1.085 
-0.07 ±0.03 -0.06 

< 0.59 0.15 


2.99l[!:J? 3.06 
1.111+°:°^? 1.069 
-0.03 ±0.07 -0.01 
-0.021 ±0.032-0.032 
< 0.63 0.22 
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here that A'^^ approximates the number of e-folds since the time when the scales associated 
with ko first crossed the horizon till the end of inflation. Thus the total number of e-folds 
since the time when the largest observable scale, kobs, crossed the horizon is ~ A'^^ + 3. 
This prior choice incorporates the theoretical prejudice that the total number of e-folds 
should be in the 50-60 range. 

The results of the Monte Carlo Markov chain (MCMC) analysis for the LOGJF and 
LOG^ cases are given in Figure [3] and are summarized in Table [2J We give 1-dimensional 
regions encompassing 68% of probability for well-determined parameters; robust upper 
bounds for parameters whose detailed constraints are parametrization-dependent; and 
best-fit values. The Table also gives posterior ranges and best-fit values for the corre- 
sponding expressions for the tilt, running and tensor-to-scalar ratio at the pivotal scale, 
no, dn/d\nk\o and ro, respectively. These have been obtained by using lowest order 
expressions in terms of the slow-roll parameters, equations (12.41) . (12.71) and (12. 8p . 

As anticipated, constraints on and lnP° are quite tight, and we have checked that 
they are almost independent of the choice of prior by performing a run with priors fiat 
in {\nN^,\nq,lnP^} instead. It is interesting to notice that in the LOGJF case the ID 
68% (95%) posterior region (2-tails) is approximately 21 < iV° < 46 (16 < iV° < 81), 
even though the mean is somewhat lower, at around 33 e-folds, and a best-fit around 26 
e-folds. This result is close to the theoretical prejudice A''^ ~ 50, required to solve the 
horizon problem. Thus assuming a fiat tree-level potential for the infiaton, the observed 
shape of the power spectrum appears to automatically point to model parameters giving 
a very sensible number of e-folds, in particular given the heavy tail of the probability 
distribution function (pdf) for large N^. This is not trivial at all: in principle any value 
for could have emerged from the analysis. 

The values of the spectral index and its running at the pivotal scale are easily derived 
from the parameters of the fit and are also given in Table [2j Note that dn/d\nk\Q is very 
small (of order ~ 10~^), as expected from the relation (I4.10p . So the LOG scenario is 
indeed close to the n = constant limit. The value of n at the pivotal scale, no, is directly 
related to the value of by equation (15. 3p . leading to the values of no quoted in the Table. 
The LOGJF best fit value, no = 0.961 (corresponding to = 25.6), coincides with the 
value obtained assuming n = constant and negligible running (Table 1). It is interesting 
to note that although the LOG^ and n = constant fits are very similar, they are not 
identical, and indeed LOGJF gives a slightly better fit, as can be checked by comparing 
the best-fit likelihood values (also compare Figure H]). Furthermore, if future CMB and 
LSS data favour a value of no closer to ~ 0.98, the value of will come out even closer 
to the theoretically preferred value, ~ 50. The upper bound on g is a consequence 
of the measured tilt and of the physical boundaries imposed on the potential parameters 
[Recall from the discussion after equation (15. lip that we expect q < {1 — n)/2.] For the 
reasons explained below, the pdf for q depends on the prior chosen, and therefore we do 
not show it in Figure [31 However the upper bound is robust with respect to a change of 
priors, and therefore we chose to report only this value. The tensor contribution remains 
negligible, below the level ~ 10"'^, since from (14.90 the value of the tensor-to-scalar ratio 
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Figure 3: ID marginalized probability distributions for the well-constrained parameters 
in the LOG scenario (compare Table [2]). Black curves are for the case with a flat prior on 
A^g , while red is for the case where a Gaussian prior around A'^^ = 50 has been enforced. 
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Table 2: Marginalized 68% regions and best — fit values for the class of models 
LOG (small-coupling regime) for quantities that are well-determined and essentially 
prior/parametrization independent. For q and the tensor-to-scalar ratio at the piv- 
otal scale, ro, we give absolute upper limits that are a consequence of the spectral 
tilt and of physical priors on the potential parameter space, Eq. fl5.5p . These bounds 
have no confidence level attached as the precise numerical value would depend on the 
prior/parametrization choice, a consequence of the PGA component being an uncon- 
strained, degenerate direction in parameter space (see text for details). 



Model 
-2 A In/: 



LOG^ 



2.1 



ID 68^ 



Best fit 



LOGJ^ 



-0.4 

ID 68% Best fit 



Gosmological parameters 



X 10^ 



r 



[Km s-^Mpc"^] 



2.30 ±0.04 2.30 

0.107 ±0.004 0.105J 

1.044 ±0.003 1.044 

0.103 ±0.026 0.107 

74.5 ± 1.5 74.9 



2.24 ±0.07 2.23 
0.107 ±0.004 1.071 
1.042 ±0.003 1.041 



0.089+^:^^' 0.089 



-0.028 
-0.030 

73.2 ± 1.9 72.9 



Power spectra parameters 



ln(P° X 10 



10^1 



3.14 ±0.05 
48.3 ±5.1 



3.14 
48.5 



3.16 



+0.05 
-0.06 



oo 4-f 12.3 



3.17 

25.6 



< 0.04 for any parametrization 



Derived power spectra parameters 



no 



dn/dlnA;|^ x 10^ 



0.979 ± 0.002 0.979 
-0.45 ±0.09 -0.43 



0.964+°:°}^ 0.961 
-1.42 ±0.94 -1.54 



< 4 X 10 ^ for any parametrization 



Potential parameters and PGA components 



(f)o/Mp 
£1 

£3 



< 4 X 10~^° for any parametrization 

< 1 X 10^^^ for any parametrization 

< 1 (from prior, Eq. (15. 5p ) 



0.00 ±0.02 
0.52 ±0.18 



0.01 
0.67 



0.00 ±0.02 
0.00 ±0.48 



0.00 
-0.20 



essentially unconstrained 
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at the pivotal scale is 

-o-4|^- (6.1) 

e 

Consequently the upper bound on q corresponds to an order of magnitude smaller bound 
on tq. 

The fact that we can extract two measured quantities in this scenario (the tilt and 
normalization) from three model parameters [either (15.11) or (15. 2p ] means that we expect 
a strongly degenerate direction in the primordial power spectrum parameter space. In 
fact, the constraints coming from the data define a region shaped as a long solid cylinder 
in the 3D subspace spanned by (15.21) . Since this cylinder is not aligned with the potential 
parameters direction, if one tries to convert limits on (15. ip into limits on the potential 
parameters (15. 2p one unavoidably picks up the degenerate direction, i.e. along the axis 
of the cylinder. This means that while in the set (15. ip the constraints on A'^^,lnP° are 
robust with respect to a change in the parametrization of the problem (since all the 
parametrization-dependency is dumped into g), it is impossible to translate these into 
parametrization-independent results for the potential parameters (15. 2p . 

However, one can still define well-constrained (and parametrization-independent) di- 
rections in the subspace spanned by (15. 2p by performing a Principal Component Analysis 
(PCA), i.e. by rotating into a new coordinate system aligned with the degenerate di- 
rection. We therefore consider the covariance matrix C in the subspace spanned by the 
reduced variables C, = (Inpo, Ing, ln0o), where hats indicate that the variables have been 
shifted by their posterior mean and normalized to their posterior standard deviation. 
Then the PCA vector e is given by 

£ = f/C, (6.2) 
where U is the 3D rotation matrix that diagonalizes C: 

CCC = e'Ae , (6.3) 

and A = diag(Ai, A2, A3) is the matrix of eigenvalues, whose square roots give the error 
along the directions defined by e. The matrix U is numerically given by 

-0.68 0.57 

U = I 0.80 0.05 -0.60 I , (6.4) 
0.73 0.56 

and = 0.02, \/\2 = 0.48 while a/As ^ 1, showing that £3 is indeed the degenerate 
direction. We have checked that the constraints on (ei, 52) are largely independent on the 
chosen parametrization. 

The 1-dimensional marginalized probability distributions for the well-constrained pa- 
rameters in the problem are shown in Figure [31 We do not show the probability distribu- 
tions for the non-primordial cosmological parameters as they are mostly very similar to 
the standard scenario, nor do we plot the pdf 's for the parameters for which we have only 
upper limits (g, tq, £3) since their distribution depends on the parametrization employed. 

Turning now to the LOG^ case, which imposes a theoretically motivated prior on the 
number of e-folds, we notice that the prior enforces ~ 50. This means that the model 
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essentially loses one further parameter, and therefore the best-fit log-likelihood is slightly 
worse (see Table [2]). In fact, the LOG^ fit has basically just one free parameter for the 
power spectrum (namely the normalization lnP°), since q is almost irrelevant. Still, it 
gives an excellent best fit to observational data! Also, enforcing 50 e-folds results in a 
very strong prediction for the tilt to be Hq ~ 0.98 (compare Figure [3] and the tightness 
of the posterior probability for no, red curve), while both the tensor contribution and the 
running are predicted to be very small. 

In Figure H] we plot the CMB temperature power spectrum for the best fit models 
discussed here along with the compilation of the data used. 

6.3 Logarithmic regime and non— renormalizable operator 
(LOG + NRO) 

In this scenario the five parameters we use to describe the primordial spectrum are 

Plog+nro = {\nPlNlq,A,N}. (6.5) 

As discussed in Section 15.21 the meaning of the first three ones is similar to those of the 
LOG case, and we impose flat priors on InPj' and q. In analogy, we consider two types 
of fits: LOG+NROJ^ (with a fiat prior on iV°) and LOG+NRO^ (with a Gaussian prior 
centered at A^^^ = 50 with standard deviation of 5). Since A is expected to be of order 
unity or less according to the discussion above, it is appropriate to use a fiat prior on A 
between and 1. 

Let us recall that A^ determines the order of the NRO, equation (13.61) . As discussed 
in Section [521 ^6 expect it to be within 1 to 0(10). One could imagine treating A^ as a 
free parameter and trying to derive a posterior bound on it from the data. However it is 
technically difficult to ensure that the MCMC is correctly performed across disjoint regions 
of the parameter space (since A^ is an integer, using it as a free parameter effectively gives 
A^ separated patches across which it is very difficult to sample). Furthermore, NRO's with 
different values of A^ are best considered as different models, since the underlying physics 
is likely to be different. Therefore distinguishing between values of A^ can be regarded 
as a model selection task, rather than a parameter constraint exercise. For this reason it 
is more instructive to consider two separate cases which are representative of the general 
behaviour at low (A^ = 2) and large (A^ = 10) values of A^. Parameter constraints from 
CMB and LSS data, which are discussed next, are summarized in Table [3] for the N = 2 
case and in Table H] for the A^ = 10 case. 

Starting from the N = 2 case, we find a strong upper bound on q, which reflects the 
theoretical considerations exposed above and is a consequence of the physically motivated 
prior (15.141) . As a consequence, the tensor contribution is always negligible. The number 
of e-folds for the LOG+NROJ^ case (A^ = 2) is = 14.5 ± 3.5 at 68%, becoming 
10.1 < < 25.8 at 95%, which is too small to solve the horizon problem. Meanwhile, the 
parameter A is rather tightly constrained, A = 0.60~^q^I. These results can be intuitively 
understood in the following way. As discussed in Subsection 4.2, the presence of the NRO 
increases dn/dlnk. This effect is maximal at low k. The lower A^ is, the more gradual 
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Figure 4: CMB temperature power spectrum for the best-fit model parameters for the 
standard parametrization, the LOG and the LOG-fNRO scenarios. The bottom panel 
shows the details of the large-scale region. 
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Table 3: As in Tabled! but for the the class of models referred in the text as LOG+NRO, 
for N = 2. Upper or lower bounds at the specified confidence level are understood to be 
1-tail limits. 



Model 
-2Aln£ 






LOG+NRO^ 


LOG±NROJ^ 


2.4 

ID 68% Best fit 


-2.7 

ID 68% Best fit 


Cosmological parameters 


fie h'^ 

r 

Hq [Km s^^Mpc^i] 


2.32 ±0.05 2.31 
0.107 ±0.004 0.107 
1.0441°:°°^ 1.044 
0.112 ±0.026 0.010 
74.7l}| 75.0 


2.18 ±0.07 2.16 
0.108 ±0.004 0.109 
1.041 ±0.003 1.042 
0.095 ±0.030 0.102 
73.3l^:^ 71.6 



Power spectra parameters 



ln(po X 10 



10^ 



e 

A 
Q 



3.14 ±0.05 3.13 
47.0 ±5.1 49.6 
0.46 ±0.11 0.27 



3.15 ±0.05 
14.5 ±3.5 



3.15 
11.7 



0.60+°:°^ 0.66 



< 5 X 10 ^ for any parametrization 



Derived power spectra parameters 



no 

dn/dlnA;|^^ x 10^ 
d^n/dln^A;!^,^ x 10^ 



0.987 ±0.007 0.981 
-0.45 
0.02 ±0.02 -0.05 



-0 11+0-11 

'^•-'--'--0.05 



1.001^^:^1^ 1.027 



0.016 
0.048 

> -6.1 (95%) -4.28 

< 14.2 (95%) 7.93 



< 3 X 10 ^ for any parametrization 



Potential parameters and PGA components 



0o/Mp 
M/Mp 

£3 
£4 



< 1 X 10"^*^ for any parametrization 

< 7 X 10"^'^ for any parametrization 

< 1 X 10~^ for any parametrization 

< 1 (from prior, Eq. (15.141) ) 



0.04 ±0.02 0.06 
0.28 ±0.13 0.05 
-2.3 ±0.4 -3.1 



0.00 ±0.02 
0.01 ±0.08 
0.01 ±0.72 



0.00 

-0.02 
0.87 



essentially unconstrained 
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is the decrease of dn/dlnk with k. In the N = 2 case the value of the running of the 
spectral index is fairly constant in the region of k accessible to observations. Thus the 
model (for not too small A, which would lead back to the LOG scenario) approximately 
resembles the dn/dlnk = constant standard parametrization. We empirically know from 
the WMAP analyses [1,2] (and our own analysis in this paper) that for this standard 
parametrization the value of n at k = ko cannot be very far from n = 1. This implies 
from Equations. (I4.19p . (I4.18P and the smallness of q that A cannot be far from A ~ (2A^+ 
3)-i/(^+2)^ which explains the value A ~ 0.6. The running dn/dlnk is then determined 
by [see equation (ICTl) ]. Not surprisingly, the preferred value for the running turns 
out to be consistent with the one from the standard parametrization (compare Tables [H 
and [3]), which corresponds to the value of quoted above. This is also consistent with 
our discussion of the (too small) number of e-folds in the standard parametrization. Some 
of these features are illustrated in Figure [5] (left panel), which shows the interplay of A'^^ 
and A and their impact on the spectral index; and Figure O which shows the best -fit 
n{k) (left panel) and the curve corresponding to the posterior mean, alongside with the 
favoured 95% posterior region of n{k) for N = 2. The corresponding Ps{k) is shown in 
the right panel. 

Note that the number of free parameters is essentially the same for both the LOG+NRO 
case and the constant running parametrization: 4 for the latter (3 if r is set to zero) and 
5 for LOG+NRO JF (among which q is almost irrelevant and has been fixed), and their 
best -fit log-likelihoods are similar. We comment further on this in Section [71 

Further enforcing a sufficient number of e-folds by imposing a Gaussian prior on A''^ 
(LOG+NRO^ case in Table [3]) results in a worsening of the quality of fit (an increase 
of minus twice the best-fit log-likelihood by 2.1 with respect to the standard power-law 
case). This is because a larger A^g for a given A implies a spectrum closer to scale- 
invariance, which is rather strongly disfavoured by data (for example, [14] reports an 
evidence of 17:1 against a scale invariant spectrum). 

As observed in the LOG case, also in the LOG+NRO case the strong degeneracy 
among the potential parameters makes it impossible to robustly translate the constraints 
on Plog+nro into prior-independent constraint for the potential parameters, (13. 5p . As we 
have done above, we can still define well constrained directions in the subspace spanned 
by the parameters ( = {In pQ,lnq,ln(j)Q,lnM) (hats indicate that the variables have been 
shifted by their posterior mean and normalized to their posterior standard deviation). 
The eigenvalues of the 3 well-constrained directions are now = 0.02, = 0.08, 
^/Aa = 0.72 while ^ 1, and the corresponding rotation matrix is 



U 



/ 0.49 


-0.56 


0.43 


0.51\ 


-0.84 


-0.05 


0.23 


0.53 


0.31 


0.79 


0.04 


0.54 


0.04 


-0.26 


-0.87 


0.42/ 



(6.6) 



Figure [7] shows marginalized 1-dimensional posterior distributions for some of the well 
constrained parameters in the LOG+NRO (A^ = 2) scenario. Since the value of the tilt, 
running and running of the running at the pivotal scale are not really representative of the 
functional form of n{k) in this case, we do not show pdf 's for those quantities (even though 
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Figure 5: Dependence of n{k) on the parameters (blue, from thin to thick A'^^ = 
10,20,30,50, for fixed A = 0.6 in the = 2 case, left, and for fixed A = 0.77 in the 
= 10 case, right) and on A (red, from thin to thick A = 0.2, 0.4, 0.6, 0.8, for fixed 
A'^ = 15 in both panels) in the LOG+NRO scenario with A^ as indicated. 



their constraints are reported for completeness in Table [3]). A more faithful representation 
is actually given in Figure [61 

Let us now turn to the A^ = 10 case. The main difference with the previous N = 2 
case is that the effect of NRO is much more pronounced on large scales. In particular 
the running dn/dlnk can now be quite large at very low k and decrease very quickly, 
converging to the LOG scenario. Thus, the scenario is qualitatively different from the 
constant running standard parametrization. As for N = 2, the value of q is small and 
below the theoretically expected bound given by equation fl5.23p . which translates again 
into a negligible tensor contribution (ro < 3 x lO"'^). The preferred value of A is still 
approximately determined by the empirical condition no ~ 1. From Equations. (14.191) . 
KT8\i . this translates into A ~ (2A^ + 3)-^/(^+2) _ o.77. The fact that the model rapidly 
converges to the LOG scenario allows to increase the number of e-folds to values similar 
to those of the LOGJF case. More precisely, in the LOG+NROJ-" case, the probability 
distribution for the number of e-folds has a heavier tail for large values of N^, and the 
2-tails posterior 68% (95%) region is given by 18.8 < < 37.1 (15.8 < < 69.5), 
which appears to be solving the horizon problem within 2a. 

Some of these features are illustrated Figure [Hi Note in particular that the running of 
running for the best fit (yellow line in Fig. [8]) is positive and sizeable. The quality of the fit 
is similar to the constant running case for the standard parametrization (though the actual 
shape of the spectrum is quite different!), although slightly worse. This is not surprising. 
From subsection 6.1 we know that data prefer a negative second derivative, which cannot 
be achieved in the LOG+NRO scenario at low k [see discussion after equation fl4.19p . 
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Table 4: As in Table [3] but for = 10 in the LOG+NRO scenario. We did not perform 
in this Principal Component Analysis as for N = 2. 



Model 
-2Aln£ 



LOG+NRO^ 



2.1 



ID 68% 



Best fit 



LOG+NROJP" 



-1.8 

ID 68% Best fit 



Cosmological parameters 



ilb X 102 
fie h'^ 

T 

Ho [Km s~^Mpc"i] 



2.31 ±0.05 


2.30 


2.23 ±0.07 


2.19 


0.107 ±0.004 


0.107 


0.107 ±0.004 


0.108 


1 n44+°-°°2 

i.U44_Q 


1.044 


1.041 ±0.003 


1.041 


0.105 ±0.026 


0.104 




0.93 


74.6 ± 1.6 


74.5 


73.1 ± 1.8 


72.1 



Power spectra parameters 



ln(po X 10^°) 

e 

A 



3.14 ±0.05 
48.2 ±5.1 



0.52 



1-0.19 
-0.16 



3.14 
48.7 
0.54 



3.15 ±0.05 



28.5 
0.57" 



+8.6 
-9.6 

-0.17 
-0.20 



3.16 
17.5 
0.77 



< 2 X 10 for any parametrization 



Derived power spectra parameters 



no 

dn/dln/i;|^ x 
d^n/dln^A;! 



10^ 
X 10^ 



0.982 



-1-0.002 
-0.004 



-0 n+0-09 

'-'•-'--'--0.06 
U'-'-'J-O.lS 



0.980 

-0.05 

-0.03 



0.971 ±0.016 1.002 
> -3.4 (95%) -4.17 
< 17.8 (95%) 23.7 



< 3 X 10 for any parametrization 



Potential parameters 



M/Mp 



< 1 X 10"^^ for any parametrization 

< 1 X 10~^^ for any parametrization 

< 0.2 for any parametrization 
< 1 (from prior, Eq. (15.141) ) 
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Figure 6: Preferred shape of the spectral index n{k) (left) and the corresponding power 
spectrum (right) from CMB and LSS data (at 95%, red curves) in the LOG+NRO scenario 
for N = 2. The yellow line shows the best-fit value, the green line the posterior mean 
while the cyan line is the best-fit further imposing a Gaussian prior on the number of 
e-folds (LOG+NRO^ scenario). The dotted blue lines represent for reference the best- 
fit power spectra in the standard parametrization with tilt only, with running and with 
running of the running (from top to bottom on the right-hand side of the n{k) panel, 
from top to bottom on the left-hand side of the Ps{k) panel, compare Table [1]). 



However, as explained there, if the spectral index really runs, a positive second derivative 
can be much more satisfactory from the physical point of view, in particular to produce 
a reasonable number of e-folds. This is precisely the case here. 

If one enforces a Gaussian prior around iV° = 50 (LOG+NRO^ case), then the best-fit 
spectrum becomes again featureless (light blue, solid line in Figure [H]) but with a smaller 
tilt than the standard parametrization (no = 0.980 for the N = 10 LOG+NRO^ case), 
which in turn means that the goodness of fit becomes worse than the standard case (see 
Table H]). Actually the model becomes in this case quite similar to the simpler LOG^ 
scenario. 

Finally, for values of between 2 and 10 we have found that the behaviour is inter- 
mediate between the cases discussed in the text. 

The above discussion shows that there are cases where the standard Taylor expansion 
of n{k) fails to capture the physics of the models. Generally, the LOG+NRO scenario 
predicts a running of n which is stronger on large scales. This can only be recovered with 
several terms in the Taylor expansion, which results in a higher number of free parameters 
in the fit. On the other hand, the functional form of n{k) in the LOG-I-NRO scenario 
implies the positiveness of the second derivative in most of the parameter space, unlike 
the standard fit. These facts make it impossible to use the results of the standard fits 
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Figure 7: ID marginalized probability distributions for some well-constrained parameters 
in the LOG+NRO scenario, for = 2 (compare Table [3]). Black curves are for the case 
with a flat prior on A'^^, while red is for the case where a Gaussian prior around A'^^ = 50. 
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Figure 8: As in Figure [6] but for the LOG+NRO scenario with = 10. 



to constrain the LOG+NRO scenario: a direct comparison of the model with the data 
becomes necessary (as we have done here). This situation could easily apply to other 
theoretical models, as well, and therefore great caution is necessary when interpreting 
generic constraints on the coefficient of the standard Taylor expansion in terms of specific 
physical models. 



7 Model comparison 

In the previous section we have presented parameter constraints for each class of model, 
namely the phenomenological model in the standard parametrization and the more phys- 
ically motivated LOG and LOG+NRO scenarios. Assessing the relative performance of 
the three models is a model comparison question, to which we now turn our attention. 

In the traditional frequentist approach to statistics, model comparison is tackled in 
terms of hypothesis testing: for example, we might ask whether the improvement in the 
best-fit likelihood in terms of the effective Ax^ = — 2Aln£ when adding a running to 
the tilt is "significant" enough to warrant the inclusion of a non-zero running. There are 
however several reasons why answering this question is far from trivial. A technical reason 
is that the usual rule of thumb of "Ax^ per extra degree of freedom" can only be applied 
if certain regularity conditions are met, and in particular only if the extra parameter 
for nested model does not lie on the boundary of the parameter space (see [35] for an 
astrophysical example and references therein). So for example, the Ax^ criterion could 
not be applied to compare the quality of fit of the LOG model with that of LOG+NRO, 
since the former is obtained from the latter by setting A = 0, and A < is not allowed. 

Another, more fundamental aspect has to do with the meaning and interpretation of 
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frequentist hypothesis testing. As discussed in detail in [13], frequentist hkehhood ratio 
tests assume the hypothesis Ti. is true and give the probabihty of observing data d as 
extreme or more extreme than what has actually been measured. This is a statement on 
the probability of the data assuming a hypothesis Ti to be true (which in Bayesian terms 
amounts to the choice of a model, Ai), i.e. frequentist hypothesis testing gives P{d\T-C). 
But this is not the quantity one is usually interested in, which is actually P{Ai\d), the 
probability of the model Ai given the observations, which can only be obtained by using 
Bayes theorem to invert the order of conditioning. For this reason, model selection is an 
inherently Bayesian question [14]. 

Bayesian model selection is based on the computation of the model likelihood P{d\Ai) = 
S{M.) (also called "evidence"), which is the normalization constant in the denominator 
of Bayes theorem (see [14,36] for details) obtained by averaging the likelihood P{d\d,Ai) 
over the prior P{6\Ai) in the parameter space 6 of the model 



From the model likelihood one obtains the model probability given the data by using 
once more Bayes theorem, P(Ai\d) oc P{A4)S{Ai), where P{Ai) is the prior probabil- 
ity assigned to the model (usually taken to be noncommittal and equal to 1/Nm if one 
considers different models). When comparing two models one usually computes the 
Bayes factor B12, given by the ratio of the evidences between the two models: 



The Bayes factor thus gives the factor by which the relative odds among two models have 
changed after the arrival of the data. As a simple calculation shows for the case of Gaussian 
likelihood and prior, equation (17. ip contains both a likelihood ratio term which rewards 
better fitting, and an "Occam's razor" term that disfavours unnecessary model complexity, 
defined in terms of useless parameters (see [37] for a discussion of model complexity). The 
"best" model is one that combines good fitting with model predictivity. Bayes factors are 
usually interpreted against the Jeffreys' scale for the strength of evidence, which we qualify 
as follows: "weak evidence" for | lni?| < 2.5, "moderate evidence" for 2.5 < | lni?| < 5.0 
and "strong evidence" for |lni?| > 5.0. The computation of the model likelihood is in 
general a numerically difficult task, as it involves a multi-dimensional integration over 
the whole of parameter space. Furthermore, a prior dependence is (correctly) built into 
the method, as the Occam's razor term depends on the ratio of the prior to posterior 
volume, which gives the amount of "wasted" parameter space of the model. Therefore it 
is problematic to evaluate the Bayes factor unless one has a physically motivated way of 
setting the prior volume. 

The difficulty of using a fully Bayesian approach to the comparison of the standard 
Taylor series parametrization with either the LOG or the LOG+NRO scenario is that 
the former represents a purely phenomenological fit to the data, while the LOG and 
LOG+NRO models are physically motivated. In particular, setting a prior on the poten- 
tial parameters of the LOG and LOG+NRO models is not comparable to setting a strictly 
phenomenological prior on the quantities of direct relevance for the fit, i.e. the spectral 
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tilt, the running, etc., in the standard parametrization. The Occam's razor effect which 
rewards highly predictive models does not work properly if we do not compare like with 
like, i.e. if we are unable to set priors on the parameter space of the phenomenological 
parametrizations used for the fit. Since the standard parametrization is by construction 
phenomenological, it cannot be directly compared using the Bayesian evidence to the 
LOG and LOG+NRO scenarios. 

However, we can still draw some interesting, partial conclusion from a Bayesian ap- 
proach. In ref. [13], a method was presented to derive upper bounds on the Bayesian 
evidence for nested models, called "Bayesian calibrated p-values" , that is useful in cases 
such as this where there is only a very loose physical basis to assign priors to phenomeno- 
logical quantities in the fit (here, the various terms in the expansion of the potential). 
This allows to assess whether extra parameters are unnecessary within the framework of 
nested models, as it gives the Bayes factor which (under mild assumptions) maximizes 
the evidence in favour of the more complex model (i.e., with more terms in the Taylor 
expansion). If this turns out to be not very strong, then one can confidently conclude 
that the extra parameters are not needed. 

Table 5: Summary of model comparison statistics. Wherever the Bayes factor is given, the 
notation In Bij indicated the Bayes factor between model i and model j (with In Bij > 
favouring model i). An overbar indicates a prior-independent upper limit obtained using 
the Bayesian calibrated p-values method. The quantity n gives the number of effective 
parameters in the model. 



Model 


Ax' 


n 


> 50 e-folds? 


Bayes factor 


Notes 


Standard parametrization 


no running 




0.0 


3(+l) 


ad hoc 






with running 




-3.4 


4(+l) 


ad hoc 


lni?2i = 0.7 


No evidence 


running of running 




-4.4 


5(+l) 


ad hoc 


\nB32 = 0.0 
lni?3i =0.4 


No evidence 
No evidence 


LOG models 


LOGJ^ 




-0.4 


~ 2 


yes (at 2a) 






LOG^ 




+2.1 


~ 1 


yes 


In 545 = 0.4 




LOG+NRO models 


LOG+NROJ^, = 


2 


-2.7 


~ 3 


no 




Excluded 


LOG+NRO^, N = 


2 


+2.4 


~ 2 


yes 




Disfavoured 


LOG+NROJ^, = 


10 


-1.8 


~ 3 


yes (at 2cr) 


\nBs4 ~ 0.0 


Au par LOGJF 


LOG+NRO^, A^ = 


10 


+2.1 


~ 2 


yes 


Infigs < -2.3 


Disfavoured 
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Table [5] summarizes some relevant model comparison statistics. Focusing first on the 
"standard parametrization" section, we have employed the method of ref. [13] to derive 
a prior-independent upper bound on the Bayesian evidence in favour of extra terms in 
the Taylor expansion. The maximum Bayesian evidence in favour of a running is only 
\nB21 = 0.7 (compared to a model with just a spectral tilt), which falls short of even 
the "weak evidence" threshold. The maximum evidence in favour of a third term in 
the Taylor expansion is even weaker. We can therefore conclude that, for the standard 
parametrization, present data do not require any higher order terms than a spectral tilt 
(for which ref. [13] found a maximum evidence of In i? = 2.9 compared to a scale-invariant 
spectrum). Notice that in this phenomenological approach the number of e-folds has to 
be added in by hand as an extra parameter of the model (although it would be derivable 
given a specific enough model for the inflaton potential), indicated as (+1) in the column 
giving an approximate value of the effective number of parameters in the model. 

Regarding the LOG class of models, the goodness of fit of the LOGJF case is similar 
to the one of the simple tilt model. Although the number of free parameters of the LOG 
scenario is 3, the parameter q is irrelevant to the fit and therefore the effective number 
of parameters is closer to 2. A more precise counting of the effective parameters could 
be achieved using the notion of Bayesian complexity [37], but this is not required in the 
context of the present discussion. It is interesting to notice that the LOG scenario also 
solves the horizon problem (within 2a of the posterior mean) with an extreme economy of 
free parameters. The LOG^ case dispenses with one further parameter (as becomes 
almost fixed to ~ 50) and the upper bound on the Bayesian evidence in favour of 
LOGJ-" indicates that the difference of Ax^ = 2.5 between the flat and Gaussian prior on 
iVg is not strongly significant. 

The LOG model can be considered nested within the LOG+NRO class of models, 
with the former formally obtained from the latter by setting A = 0. For N = 2, the 
LOG+NROJF model falls short of achieving the necessary number of e-folds, and for this 
reason it must be excluded, even though its quality of fit is comparable to the standard 
case with constant running. The LOG+NRO^ case has one extra parameter (for fixed 
A^) than LOG^, and a best-fit value which is actually slightly worse, a consequence of 
the Gaussian prior forcing the posterior distribution around a value of which is not 
strongly favoured by the data. Hence we can conclude that LOG+NRO^ (A^ = 2) is 
disfavoured with respect to LOG^ and LOGJF since it is unable to achieve a better fit 
even with one extra parameter. 

The LOG+NROJF with A^ = 10 has a better fit than the LOG scenario and it achieves 
a sufficient number of e-folds within 2a. The method of Bayesian calibrated p-values 
cannot be used to compare the two models because the LOG model (obtained by setting 
y4 = in the LOG+NRO model) lies at the boundary of parameter space. However, we 
can still roughly estimate the Bayes factor between LOGJF and LOG+NROJ-" by taking 
a prior width on the extra parameter A of order unity (as motivated by the theoretical 
expectations presented in Section [5l2|) and using that, for nested models, the Bayes factor 
in favour of the simpler model is approximately (see equation (9) in [14]) 

ln5~/-AV2, (7.3) 
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where / is the logarithm of the ratio of the prior to posterior volume (the information gain) 
for the extra parameter and A is the number of sigmas discrepancy between the likelihood 
peak and the value of the extra parameter under the nested model (here, A = 0). Using 
the values in Table H] one obtains A ~ 2.8 and I ~ 1.6 and thus InS ~ —2.3, which 
would weakly favour the LOG+NROJF model. However, one has to bear in mind that 
the parameter has been fixed to a value picked among a range of order 10 possible 
values — hence one has to factor in an extra Occam's razor effect coming from the fact 
that = 10 is one of about 10 possible choices for A^. Hence Ini? has to be increased 
by about a factor In 10 = 2.3, which brings the final odds between LOG+NROJF and 
LOGJF to unity (i.e., InB ~ 0). Finally, the LOG+NRO^ case has the same quality of 
fit of the LOG^ case and one extra parameter. The Occam's razor term from the choice 
of alone would disfavour LOG+NRO^ by a factor InB = 2.3 with respect to LOG^, 
so even without computing the precise Bayes factor we can conclude that this scenario is 
disfavoured. 

In conclusion, a model comparison approach singles out the LOG scenario and the 
LOG+NROJF (A^ = 10) model as the most viable cases in light of the present data. 
This kind of considerations could be extended to compare this class of models with other 
inflationary scenarios, once they have been suitably parametrized in terms of fundamental 
variables. However, a direct comparison with a phenomenological approach such as the 
standard Taylor expansion of the spectrum is not feasible due to the lack of predictivity 
of the latter. The Bayesian evidence still concludes that no higher-order term than the 
tilt is presently required in the series. 

Finally, we emphasize that the LOG and LOG+NRO models predict tensor contribu- 
tions that are generally very small and will be largely undetectable. The most optimistic 
case is the LOG, where the upper bound is of order tq ~ 10~^, which might be just within 
reach of future B-modes observations. Conversely, a detection of tensor modes above 
~ 10^^ would disprove the scenario of fiat tree-level inflationary potentials. 

8 Conclusions 

In this paper we have compared a broad and physically well-motivated class of inflationary 
models with CMB and LSS observational data. Namely, we have considered models with 
flat tree-level potentials, which typically appear in supersymmetric theories, where Vj.^^^"^ 
ordinarily has plenty of accidental flat directions. These models, beside being very well 
motivated from the physical point of view (on a similar footing to monomial potentials), 
lead to very model-independent cosmological predictions. The reason is that the potential 
derivatives V, V",... arise from the radiative corrections to V, which has a characteristic 
logarithmic dependence on the inflaton field. This scenario has been labelled "LOG" 
throughout the paper. In addition, we have considered the possible presence of new 
physics beyond a certain high-energy cut-off. This physics does not need to respect the 
flat directions of the "low-energy" theory, and thus it will show up as non-renormalizable 
operators (NRO) in the inflaton field, which will be dominated by the lowest-order one. 
This modified scenario (labelled "LOG-I-NRO") is also very well motivated and still quite 
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model-independent. 

We have studied the performance of these scenarios when compared with CMB and 
LSS. We have made first a detailed study of the features of these models, working out 
both numerically accurate results and approximate analytical expressions for Ps{k), Pt{k) 
and other relevant quantities, such as the spectral index, n{k), as a function of suitably 
defined model parameters. We also discussed the number of independent parameters and 
the theoretical and phenomenological constraints on them (to be imposed a priori). As 
a matter of fact, one (combination) of the parameters is almost irrelevant, which makes 
these models even more predictive. Another parameter is essentially the number of e-folds, 
N^, since the time when the largest observable scales crossed out the horizon until the end 
of inflation. This allows to perform the fits in a twofold way: either leaving A^^^ free, and 
let the data determine its value, or imposing a prior on its value according to the usual 
theoretical prejudice {N^ = 50 — 60). Both approaches (labelled and Q respectively) 
are interesting and complementary. 

In the analysis we also study the performance of standard parametrizations of the 
power spectrum, based on a Taylor expansion of lnP5(A;) and In Pt{k) around an (arbi- 
trary) pivotal scale, Iu/cq. At first (second) order these parametrizations correspond to 
a constant (constantly running) spectral index, n{k) = constant {dn/d\nk = constant). 
They have been used in reference analyses, in particular by the WMAP collaboration. 
It is important to keep in mind that, although useful, the standard parametrizations are 
not inspired by any particularly well-motivated inflationary physics. E.g. the results of 
the fit with the dn/d\nk = constant assumption are not consistent with a number of 
e-folds in the required range. Still we have also studied them (going one order beyond 
the WMAP analysis) to facilitate the discussion of the performance of the LOG and 
LOG-I-NRO scenarios. As a general comment, care must be taken to test inflationary 
models which predict a non-negligible scale dependence. In many cases the standard Tay- 
lor series parametrization of equations 02.91) and fl2.10|] cannot be accurately used in such 
a situation unless a high number of terms is taken in the expansion. 

Our main results are the following: 

• Both the LOG and LOG+NRO scenarios predict small tensor perturbations: tq < 
0(10-3). 

• The LOG scenario has essentially two parameters, Ps{ko) and N^, and implies a 
nearly constant n{k). 

• Leaving N^, free parameter (LOGJT fit), one gets 24 < < 49 (16 < A^^ < 84) 
at 68% (95%) c.l. while the corresponding spectral index is close to Uq = 0.96. This 
result is consistent with the theoretical prejudice ~ 50 — 60, required to solve the 
horizon problem, which is remarkable. In the LOG^ fit (i.e. imposing = 50 — 60) 
one gets no ^ 0.98. Note that this fit has only one parameter and still works very 
well. 

• The LOG-I-NRO scenario has two more parameters, which arise from the order and 
the suppression scale of the NRO. We have fixed the order of the NRO (denoted 4 + 
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27V throughout the paper) to two representative values (A^ = 2, 10) that reasonably 
encompass the sensible physical range. Thus in practice we are playing with just one 
additional parameter. This scenario can produce a sizeable running of the spectral 
index, and still be consistent with the data and a reasonable number of e-folds, 
especially if N is not very small. The impact of the NRO (driving a running of 
the spectral index) is relevant at small k, corresponding to the first stages of the 
inflationary period, and then it quickly converges (especially for not too small N) 
to the LOG scenario. 

• The model comparison is delicate for several technical and fundamental reasons. In 
the paper we give a fully Bayesian discussion of the relative quality of the various 
scenarios considered. Qualitatively, it can be said that the goodness of the LOG 
(NRO+LOGJF) fits is similar to the standard n = constant {dn/dlnk = constant) 
parametrization. The improvement in the goodness of the fit obtained by the inclu- 
sion of an extra parameter (as the LOG+NRO scenario implies) is not enough (with 
the present data) to justify such modification, but still it remains an interesting the- 
oretical possibility. On the other hand, a rigorous comparison between the evidences 
for LOG, LOG+NRO scenarios and for the standard parametrizations (which are 
phenomenological by construction) is not feasible. However, the prior-independent 
method of the Bayesian calibrated p-values still indicates that no higher order term 
than a tilt is required in the standard Taylor expansion. 

As a final conclusion, the LOG and LOG+NRO scenarios analyzed in this paper (based 
on flat tree-level potentials without or with the presence of extra physics) are not only very 
well motivated from the physical point of view, but they also fit remarkably well the CMB 
and LSS data, with very few parameters (the predictions are quite model independent). 
In addition they are naturally consistent with a reasonable number of e-folds. Therefore, 
they can be considered as a standard physical class of inflationary models, on a similar 
footing as monomial potentials. 
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A Some complementary formulas 



In this appendix we provide some important expressions for the case of the potential with 
NRO. Applying the following two identities of hypergeometric functions: 



2Fi(a,fe; c; z) 
3^2(0, 6, c; a + 1, 6 + 1; z) 



2Fi{b,a;c;z 
1 



b — a 

one can see, integrating equation fl2.5p . that 

1 /2 1 
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(Av;)-^-'-(iV + l)2Fi 
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^=1 
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(A.l) 



being the slow-roll approximation the only reason for the symbol of approximate equality. 
The expression f l4.15p that we use in the fits is directly obtained from flA.ip neglecting the 
irrelevant addends, i.e. all but the one that is not proportional to the small parameter q. 

The first order slow-roll parameters in this scenario are: 



1 

2qN^ $ 



q ' 2M 
(2M- 1)(A$) 



ln<I> 

- 1 



i + ^(A$r + iln$ 



where we have defined 



(A.2) 
(A.3) 

(A.4) 



to simplify the notation. Notice that in the limit of the NRO going to zero we recover 
the formulas (14. ip . 

The relation between Plog+nro and the physical parameters of the potential are given 

by: 
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